Here, I’m trying to get a better sense of how well the SODA temp data represents the SEUS. We found a significant decrease in temperature during the study period, which is not what people who are familiar with conditions in this area say has happened. Most likely, the SODA model doesn’t perform well in shallow near shore areas where survey is performed in southeast US. Therefore, we compared performance to OISST temp. Overall pattern is similar, and not representative of what’s actually happening in this survey area. Therefore, we decided to go with in situ buoy temperature values.
library(ncdf4)
library(raster)
Loading required package: sp
Attaching package: ‘raster’
The following object is masked from ‘package:dplyr’:
select
The following object is masked from ‘package:tidyr’:
extract
The following object is masked from ‘package:lme4’:
getData
The following object is masked from ‘package:data.table’:
shift
library(data.table)
Southeast US (1990–2014, n = 25)
Years downloaded one by one on May 15, 2020 https://psl.noaa.gov/cgi-bin/DataAccess.pl?DB_dataset=NOAA+High-resolution+Blended+Analysis&DB_variable=Sea+Surface+Temperature&DB_statistic=Mean&DB_tid=86107&DB_did=132&DB_vid=2423
#all file names
filenames <- list.files(here::here("Data","Temperature","OISST_NOAA_TEMP_NETCDF"), full.names = T)
OISST.sst <- lapply(filenames, brick)
OISST.sst
[[1]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1990.nc
names : X1990.01.01, X1990.01.02, X1990.01.03, X1990.01.04, X1990.01.05, X1990.01.06, X1990.01.07, X1990.01.08, X1990.01.09, X1990.01.10, X1990.01.11, X1990.01.12, X1990.01.13, X1990.01.14, X1990.01.15, ...
Date : 1990-01-01, 1990-12-31 (min, max)
varname : sst
[[2]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1991.nc
names : X1991.01.01, X1991.01.02, X1991.01.03, X1991.01.04, X1991.01.05, X1991.01.06, X1991.01.07, X1991.01.08, X1991.01.09, X1991.01.10, X1991.01.11, X1991.01.12, X1991.01.13, X1991.01.14, X1991.01.15, ...
Date : 1991-01-01, 1991-12-31 (min, max)
varname : sst
[[3]]
class : RasterBrick
dimensions : 32, 28, 896, 366 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1992.nc
names : X1992.01.01, X1992.01.02, X1992.01.03, X1992.01.04, X1992.01.05, X1992.01.06, X1992.01.07, X1992.01.08, X1992.01.09, X1992.01.10, X1992.01.11, X1992.01.12, X1992.01.13, X1992.01.14, X1992.01.15, ...
Date : 1992-01-01, 1992-12-31 (min, max)
varname : sst
[[4]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1993.nc
names : X1993.01.01, X1993.01.02, X1993.01.03, X1993.01.04, X1993.01.05, X1993.01.06, X1993.01.07, X1993.01.08, X1993.01.09, X1993.01.10, X1993.01.11, X1993.01.12, X1993.01.13, X1993.01.14, X1993.01.15, ...
Date : 1993-01-01, 1993-12-31 (min, max)
varname : sst
[[5]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1994.nc
names : X1994.01.01, X1994.01.02, X1994.01.03, X1994.01.04, X1994.01.05, X1994.01.06, X1994.01.07, X1994.01.08, X1994.01.09, X1994.01.10, X1994.01.11, X1994.01.12, X1994.01.13, X1994.01.14, X1994.01.15, ...
Date : 1994-01-01, 1994-12-31 (min, max)
varname : sst
[[6]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1995.nc
names : X1995.01.01, X1995.01.02, X1995.01.03, X1995.01.04, X1995.01.05, X1995.01.06, X1995.01.07, X1995.01.08, X1995.01.09, X1995.01.10, X1995.01.11, X1995.01.12, X1995.01.13, X1995.01.14, X1995.01.15, ...
Date : 1995-01-01, 1995-12-31 (min, max)
varname : sst
[[7]]
class : RasterBrick
dimensions : 32, 28, 896, 366 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1996.nc
names : X1996.01.01, X1996.01.02, X1996.01.03, X1996.01.04, X1996.01.05, X1996.01.06, X1996.01.07, X1996.01.08, X1996.01.09, X1996.01.10, X1996.01.11, X1996.01.12, X1996.01.13, X1996.01.14, X1996.01.15, ...
Date : 1996-01-01, 1996-12-31 (min, max)
varname : sst
[[8]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1997.nc
names : X1997.01.01, X1997.01.02, X1997.01.03, X1997.01.04, X1997.01.05, X1997.01.06, X1997.01.07, X1997.01.08, X1997.01.09, X1997.01.10, X1997.01.11, X1997.01.12, X1997.01.13, X1997.01.14, X1997.01.15, ...
Date : 1997-01-01, 1997-12-31 (min, max)
varname : sst
[[9]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1998.nc
names : X1998.01.01, X1998.01.02, X1998.01.03, X1998.01.04, X1998.01.05, X1998.01.06, X1998.01.07, X1998.01.08, X1998.01.09, X1998.01.10, X1998.01.11, X1998.01.12, X1998.01.13, X1998.01.14, X1998.01.15, ...
Date : 1998-01-01, 1998-12-31 (min, max)
varname : sst
[[10]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/1999.nc
names : X1999.01.01, X1999.01.02, X1999.01.03, X1999.01.04, X1999.01.05, X1999.01.06, X1999.01.07, X1999.01.08, X1999.01.09, X1999.01.10, X1999.01.11, X1999.01.12, X1999.01.13, X1999.01.14, X1999.01.15, ...
Date : 1999-01-01, 1999-12-31 (min, max)
varname : sst
[[11]]
class : RasterBrick
dimensions : 32, 28, 896, 366 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2000.nc
names : X2000.01.01, X2000.01.02, X2000.01.03, X2000.01.04, X2000.01.05, X2000.01.06, X2000.01.07, X2000.01.08, X2000.01.09, X2000.01.10, X2000.01.11, X2000.01.12, X2000.01.13, X2000.01.14, X2000.01.15, ...
Date : 2000-01-01, 2000-12-31 (min, max)
varname : sst
[[12]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2001.nc
names : X2001.01.01, X2001.01.02, X2001.01.03, X2001.01.04, X2001.01.05, X2001.01.06, X2001.01.07, X2001.01.08, X2001.01.09, X2001.01.10, X2001.01.11, X2001.01.12, X2001.01.13, X2001.01.14, X2001.01.15, ...
Date : 2001-01-01, 2001-12-31 (min, max)
varname : sst
[[13]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2002.nc
names : X2002.01.01, X2002.01.02, X2002.01.03, X2002.01.04, X2002.01.05, X2002.01.06, X2002.01.07, X2002.01.08, X2002.01.09, X2002.01.10, X2002.01.11, X2002.01.12, X2002.01.13, X2002.01.14, X2002.01.15, ...
Date : 2002-01-01, 2002-12-31 (min, max)
varname : sst
[[14]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2003.nc
names : X2003.01.01, X2003.01.02, X2003.01.03, X2003.01.04, X2003.01.05, X2003.01.06, X2003.01.07, X2003.01.08, X2003.01.09, X2003.01.10, X2003.01.11, X2003.01.12, X2003.01.13, X2003.01.14, X2003.01.15, ...
Date : 2003-01-01, 2003-12-31 (min, max)
varname : sst
[[15]]
class : RasterBrick
dimensions : 32, 28, 896, 366 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2004.nc
names : X2004.01.01, X2004.01.02, X2004.01.03, X2004.01.04, X2004.01.05, X2004.01.06, X2004.01.07, X2004.01.08, X2004.01.09, X2004.01.10, X2004.01.11, X2004.01.12, X2004.01.13, X2004.01.14, X2004.01.15, ...
Date : 2004-01-01, 2004-12-31 (min, max)
varname : sst
[[16]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2005.nc
names : X2005.01.01, X2005.01.02, X2005.01.03, X2005.01.04, X2005.01.05, X2005.01.06, X2005.01.07, X2005.01.08, X2005.01.09, X2005.01.10, X2005.01.11, X2005.01.12, X2005.01.13, X2005.01.14, X2005.01.15, ...
Date : 2005-01-01, 2005-12-31 (min, max)
varname : sst
[[17]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2006.nc
names : X2006.01.01, X2006.01.02, X2006.01.03, X2006.01.04, X2006.01.05, X2006.01.06, X2006.01.07, X2006.01.08, X2006.01.09, X2006.01.10, X2006.01.11, X2006.01.12, X2006.01.13, X2006.01.14, X2006.01.15, ...
Date : 2006-01-01, 2006-12-31 (min, max)
varname : sst
[[18]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2007.nc
names : X2007.01.01, X2007.01.02, X2007.01.03, X2007.01.04, X2007.01.05, X2007.01.06, X2007.01.07, X2007.01.08, X2007.01.09, X2007.01.10, X2007.01.11, X2007.01.12, X2007.01.13, X2007.01.14, X2007.01.15, ...
Date : 2007-01-01, 2007-12-31 (min, max)
varname : sst
[[19]]
class : RasterBrick
dimensions : 32, 28, 896, 366 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2008.nc
names : X2008.01.01, X2008.01.02, X2008.01.03, X2008.01.04, X2008.01.05, X2008.01.06, X2008.01.07, X2008.01.08, X2008.01.09, X2008.01.10, X2008.01.11, X2008.01.12, X2008.01.13, X2008.01.14, X2008.01.15, ...
Date : 2008-01-01, 2008-12-31 (min, max)
varname : sst
[[20]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2009.nc
names : X2009.01.01, X2009.01.02, X2009.01.03, X2009.01.04, X2009.01.05, X2009.01.06, X2009.01.07, X2009.01.08, X2009.01.09, X2009.01.10, X2009.01.11, X2009.01.12, X2009.01.13, X2009.01.14, X2009.01.15, ...
Date : 2009-01-01, 2009-12-31 (min, max)
varname : sst
[[21]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2010.nc
names : X2010.01.01, X2010.01.02, X2010.01.03, X2010.01.04, X2010.01.05, X2010.01.06, X2010.01.07, X2010.01.08, X2010.01.09, X2010.01.10, X2010.01.11, X2010.01.12, X2010.01.13, X2010.01.14, X2010.01.15, ...
Date : 2010-01-01, 2010-12-31 (min, max)
varname : sst
[[22]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2011.nc
names : X2011.01.01, X2011.01.02, X2011.01.03, X2011.01.04, X2011.01.05, X2011.01.06, X2011.01.07, X2011.01.08, X2011.01.09, X2011.01.10, X2011.01.11, X2011.01.12, X2011.01.13, X2011.01.14, X2011.01.15, ...
Date : 2011-01-01, 2011-12-31 (min, max)
varname : sst
[[23]]
class : RasterBrick
dimensions : 32, 28, 896, 366 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2012.nc
names : X2012.01.01, X2012.01.02, X2012.01.03, X2012.01.04, X2012.01.05, X2012.01.06, X2012.01.07, X2012.01.08, X2012.01.09, X2012.01.10, X2012.01.11, X2012.01.12, X2012.01.13, X2012.01.14, X2012.01.15, ...
Date : 2012-01-01, 2012-12-31 (min, max)
varname : sst
[[24]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2013.nc
names : X2013.01.01, X2013.01.02, X2013.01.03, X2013.01.04, X2013.01.05, X2013.01.06, X2013.01.07, X2013.01.08, X2013.01.09, X2013.01.10, X2013.01.11, X2013.01.12, X2013.01.13, X2013.01.14, X2013.01.15, ...
Date : 2013-01-01, 2013-12-31 (min, max)
varname : sst
[[25]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2014.nc
names : X2014.01.01, X2014.01.02, X2014.01.03, X2014.01.04, X2014.01.05, X2014.01.06, X2014.01.07, X2014.01.08, X2014.01.09, X2014.01.10, X2014.01.11, X2014.01.12, X2014.01.13, X2014.01.14, X2014.01.15, ...
Date : 2014-01-01, 2014-12-31 (min, max)
varname : sst
[[26]]
class : RasterBrick
dimensions : 32, 28, 896, 365 (nrow, ncol, ncell, nlayers)
resolution : 0.25, 0.25 (x, y)
extent : 278, 285, 28, 36 (xmin, xmax, ymin, ymax)
crs : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0
source : /Users/zoekitchel/Documents/grad school/Rutgers/Repositories/Col_Ext/Data/Temperature/OISST_NOAA_TEMP_NETCDF/2015.nc
names : X2015.01.01, X2015.01.02, X2015.01.03, X2015.01.04, X2015.01.05, X2015.01.06, X2015.01.07, X2015.01.08, X2015.01.09, X2015.01.10, X2015.01.11, X2015.01.12, X2015.01.13, X2015.01.14, X2015.01.15, ...
Date : 2015-01-01, 2015-12-31 (min, max)
varname : sst
#currently, this is in 0-360, needs to be -180-180, so I have to subtract 360Ëš
OISST.sst_adj <- lapply(OISST.sst, rotate)
this does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this functionthis does not look like an appropriate object for this function
Extract specific GPS values for SE US sampled during during study period
load(here::here("Data","Spp_observations","trawl_coordinates_1996-2000.RData"))
seus_points <- gpspointsonly_quarterdegreeround[region == "SCDNR_SEUS"]
coordinates(seus_points) = ~ lon + lat
temp_values_1990 <- as.data.table(extract(OISST.sst_adj[[1]], seus_points))
temp_values_1991 <- as.data.table(extract(OISST.sst_adj[[2]], seus_points))
temp_values_1992 <- as.data.table(extract(OISST.sst_adj[[3]], seus_points))
temp_values_1993 <- as.data.table(extract(OISST.sst_adj[[4]], seus_points))
temp_values_1994 <- as.data.table(extract(OISST.sst_adj[[5]], seus_points))
temp_values_1995 <- as.data.table(extract(OISST.sst_adj[[6]], seus_points))
temp_values_1996 <- as.data.table(extract(OISST.sst_adj[[7]], seus_points))
temp_values_1997 <- as.data.table(extract(OISST.sst_adj[[8]], seus_points))
temp_values_1998 <- as.data.table(extract(OISST.sst_adj[[9]], seus_points))
temp_values_1999 <- as.data.table(extract(OISST.sst_adj[[10]], seus_points))
temp_values_2000 <- as.data.table(extract(OISST.sst_adj[[11]], seus_points))
temp_values_2001 <- as.data.table(extract(OISST.sst_adj[[12]], seus_points))
temp_values_2002 <- as.data.table(extract(OISST.sst_adj[[13]], seus_points))
temp_values_2003 <- as.data.table(extract(OISST.sst_adj[[14]], seus_points))
temp_values_2004 <- as.data.table(extract(OISST.sst_adj[[15]], seus_points))
temp_values_2005 <- as.data.table(extract(OISST.sst_adj[[16]], seus_points))
temp_values_2006 <- as.data.table(extract(OISST.sst_adj[[17]], seus_points))
temp_values_2007 <- as.data.table(extract(OISST.sst_adj[[18]], seus_points))
temp_values_2008 <- as.data.table(extract(OISST.sst_adj[[19]], seus_points))
temp_values_2009 <- as.data.table(extract(OISST.sst_adj[[20]], seus_points))
temp_values_2010 <- as.data.table(extract(OISST.sst_adj[[21]], seus_points))
temp_values_2011 <- as.data.table(extract(OISST.sst_adj[[22]], seus_points))
temp_values_2012 <- as.data.table(extract(OISST.sst_adj[[23]], seus_points))
temp_values_2013 <- as.data.table(extract(OISST.sst_adj[[24]], seus_points))
temp_values_2014 <- as.data.table(extract(OISST.sst_adj[[25]], seus_points))
temp_values_2015 <- as.data.table(extract(OISST.sst_adj[[26]], seus_points))
plot(OISST.sst_adj[[1]][[1]])
points(seus_points)
object_names <- list(temp_values_1990, temp_values_1991, temp_values_1992, temp_values_1993, temp_values_1994, temp_values_1995, temp_values_1996, temp_values_1997, temp_values_1998, temp_values_1999, temp_values_2000, temp_values_2001, temp_values_2002, temp_values_2003, temp_values_2004, temp_values_2005, temp_values_2006, temp_values_2007, temp_values_2008, temp_values_2009, temp_values_2010, temp_values_2011, temp_values_2012, temp_values_2013, temp_values_2014, temp_values_2015)
allvalues_wide <- do.call("cbind", object_names)
allvalues_wide_latlon <- cbind(seus_points@coords, allvalues_wide)
colnames(allvalues_wide_latlon)
[1] "lon" "lat" "X1990.01.01" "X1990.01.02" "X1990.01.03" "X1990.01.04" "X1990.01.05"
[8] "X1990.01.06" "X1990.01.07" "X1990.01.08" "X1990.01.09" "X1990.01.10" "X1990.01.11" "X1990.01.12"
[15] "X1990.01.13" "X1990.01.14" "X1990.01.15" "X1990.01.16" "X1990.01.17" "X1990.01.18" "X1990.01.19"
[22] "X1990.01.20" "X1990.01.21" "X1990.01.22" "X1990.01.23" "X1990.01.24" "X1990.01.25" "X1990.01.26"
[29] "X1990.01.27" "X1990.01.28" "X1990.01.29" "X1990.01.30" "X1990.01.31" "X1990.02.01" "X1990.02.02"
[36] "X1990.02.03" "X1990.02.04" "X1990.02.05" "X1990.02.06" "X1990.02.07" "X1990.02.08" "X1990.02.09"
[43] "X1990.02.10" "X1990.02.11" "X1990.02.12" "X1990.02.13" "X1990.02.14" "X1990.02.15" "X1990.02.16"
[50] "X1990.02.17" "X1990.02.18" "X1990.02.19" "X1990.02.20" "X1990.02.21" "X1990.02.22" "X1990.02.23"
[57] "X1990.02.24" "X1990.02.25" "X1990.02.26" "X1990.02.27" "X1990.02.28" "X1990.03.01" "X1990.03.02"
[64] "X1990.03.03" "X1990.03.04" "X1990.03.05" "X1990.03.06" "X1990.03.07" "X1990.03.08" "X1990.03.09"
[71] "X1990.03.10" "X1990.03.11" "X1990.03.12" "X1990.03.13" "X1990.03.14" "X1990.03.15" "X1990.03.16"
[78] "X1990.03.17" "X1990.03.18" "X1990.03.19" "X1990.03.20" "X1990.03.21" "X1990.03.22" "X1990.03.23"
[85] "X1990.03.24" "X1990.03.25" "X1990.03.26" "X1990.03.27" "X1990.03.28" "X1990.03.29" "X1990.03.30"
[92] "X1990.03.31" "X1990.04.01" "X1990.04.02" "X1990.04.03" "X1990.04.04" "X1990.04.05" "X1990.04.06"
[99] "X1990.04.07" "X1990.04.08" "X1990.04.09" "X1990.04.10" "X1990.04.11" "X1990.04.12" "X1990.04.13"
[106] "X1990.04.14" "X1990.04.15" "X1990.04.16" "X1990.04.17" "X1990.04.18" "X1990.04.19" "X1990.04.20"
[113] "X1990.04.21" "X1990.04.22" "X1990.04.23" "X1990.04.24" "X1990.04.25" "X1990.04.26" "X1990.04.27"
[120] "X1990.04.28" "X1990.04.29" "X1990.04.30" "X1990.05.01" "X1990.05.02" "X1990.05.03" "X1990.05.04"
[127] "X1990.05.05" "X1990.05.06" "X1990.05.07" "X1990.05.08" "X1990.05.09" "X1990.05.10" "X1990.05.11"
[134] "X1990.05.12" "X1990.05.13" "X1990.05.14" "X1990.05.15" "X1990.05.16" "X1990.05.17" "X1990.05.18"
[141] "X1990.05.19" "X1990.05.20" "X1990.05.21" "X1990.05.22" "X1990.05.23" "X1990.05.24" "X1990.05.25"
[148] "X1990.05.26" "X1990.05.27" "X1990.05.28" "X1990.05.29" "X1990.05.30" "X1990.05.31" "X1990.06.01"
[155] "X1990.06.02" "X1990.06.03" "X1990.06.04" "X1990.06.05" "X1990.06.06" "X1990.06.07" "X1990.06.08"
[162] "X1990.06.09" "X1990.06.10" "X1990.06.11" "X1990.06.12" "X1990.06.13" "X1990.06.14" "X1990.06.15"
[169] "X1990.06.16" "X1990.06.17" "X1990.06.18" "X1990.06.19" "X1990.06.20" "X1990.06.21" "X1990.06.22"
[176] "X1990.06.23" "X1990.06.24" "X1990.06.25" "X1990.06.26" "X1990.06.27" "X1990.06.28" "X1990.06.29"
[183] "X1990.06.30" "X1990.07.01" "X1990.07.02" "X1990.07.03" "X1990.07.04" "X1990.07.05" "X1990.07.06"
[190] "X1990.07.07" "X1990.07.08" "X1990.07.09" "X1990.07.10" "X1990.07.11" "X1990.07.12" "X1990.07.13"
[197] "X1990.07.14" "X1990.07.15" "X1990.07.16" "X1990.07.17" "X1990.07.18" "X1990.07.19" "X1990.07.20"
[204] "X1990.07.21" "X1990.07.22" "X1990.07.23" "X1990.07.24" "X1990.07.25" "X1990.07.26" "X1990.07.27"
[211] "X1990.07.28" "X1990.07.29" "X1990.07.30" "X1990.07.31" "X1990.08.01" "X1990.08.02" "X1990.08.03"
[218] "X1990.08.04" "X1990.08.05" "X1990.08.06" "X1990.08.07" "X1990.08.08" "X1990.08.09" "X1990.08.10"
[225] "X1990.08.11" "X1990.08.12" "X1990.08.13" "X1990.08.14" "X1990.08.15" "X1990.08.16" "X1990.08.17"
[232] "X1990.08.18" "X1990.08.19" "X1990.08.20" "X1990.08.21" "X1990.08.22" "X1990.08.23" "X1990.08.24"
[239] "X1990.08.25" "X1990.08.26" "X1990.08.27" "X1990.08.28" "X1990.08.29" "X1990.08.30" "X1990.08.31"
[246] "X1990.09.01" "X1990.09.02" "X1990.09.03" "X1990.09.04" "X1990.09.05" "X1990.09.06" "X1990.09.07"
[253] "X1990.09.08" "X1990.09.09" "X1990.09.10" "X1990.09.11" "X1990.09.12" "X1990.09.13" "X1990.09.14"
[260] "X1990.09.15" "X1990.09.16" "X1990.09.17" "X1990.09.18" "X1990.09.19" "X1990.09.20" "X1990.09.21"
[267] "X1990.09.22" "X1990.09.23" "X1990.09.24" "X1990.09.25" "X1990.09.26" "X1990.09.27" "X1990.09.28"
[274] "X1990.09.29" "X1990.09.30" "X1990.10.01" "X1990.10.02" "X1990.10.03" "X1990.10.04" "X1990.10.05"
[281] "X1990.10.06" "X1990.10.07" "X1990.10.08" "X1990.10.09" "X1990.10.10" "X1990.10.11" "X1990.10.12"
[288] "X1990.10.13" "X1990.10.14" "X1990.10.15" "X1990.10.16" "X1990.10.17" "X1990.10.18" "X1990.10.19"
[295] "X1990.10.20" "X1990.10.21" "X1990.10.22" "X1990.10.23" "X1990.10.24" "X1990.10.25" "X1990.10.26"
[302] "X1990.10.27" "X1990.10.28" "X1990.10.29" "X1990.10.30" "X1990.10.31" "X1990.11.01" "X1990.11.02"
[309] "X1990.11.03" "X1990.11.04" "X1990.11.05" "X1990.11.06" "X1990.11.07" "X1990.11.08" "X1990.11.09"
[316] "X1990.11.10" "X1990.11.11" "X1990.11.12" "X1990.11.13" "X1990.11.14" "X1990.11.15" "X1990.11.16"
[323] "X1990.11.17" "X1990.11.18" "X1990.11.19" "X1990.11.20" "X1990.11.21" "X1990.11.22" "X1990.11.23"
[330] "X1990.11.24" "X1990.11.25" "X1990.11.26" "X1990.11.27" "X1990.11.28" "X1990.11.29" "X1990.11.30"
[337] "X1990.12.01" "X1990.12.02" "X1990.12.03" "X1990.12.04" "X1990.12.05" "X1990.12.06" "X1990.12.07"
[344] "X1990.12.08" "X1990.12.09" "X1990.12.10" "X1990.12.11" "X1990.12.12" "X1990.12.13" "X1990.12.14"
[351] "X1990.12.15" "X1990.12.16" "X1990.12.17" "X1990.12.18" "X1990.12.19" "X1990.12.20" "X1990.12.21"
[358] "X1990.12.22" "X1990.12.23" "X1990.12.24" "X1990.12.25" "X1990.12.26" "X1990.12.27" "X1990.12.28"
[365] "X1990.12.29" "X1990.12.30" "X1990.12.31" "X1991.01.01" "X1991.01.02" "X1991.01.03" "X1991.01.04"
[372] "X1991.01.05" "X1991.01.06" "X1991.01.07" "X1991.01.08" "X1991.01.09" "X1991.01.10" "X1991.01.11"
[379] "X1991.01.12" "X1991.01.13" "X1991.01.14" "X1991.01.15" "X1991.01.16" "X1991.01.17" "X1991.01.18"
[386] "X1991.01.19" "X1991.01.20" "X1991.01.21" "X1991.01.22" "X1991.01.23" "X1991.01.24" "X1991.01.25"
[393] "X1991.01.26" "X1991.01.27" "X1991.01.28" "X1991.01.29" "X1991.01.30" "X1991.01.31" "X1991.02.01"
[400] "X1991.02.02" "X1991.02.03" "X1991.02.04" "X1991.02.05" "X1991.02.06" "X1991.02.07" "X1991.02.08"
[407] "X1991.02.09" "X1991.02.10" "X1991.02.11" "X1991.02.12" "X1991.02.13" "X1991.02.14" "X1991.02.15"
[414] "X1991.02.16" "X1991.02.17" "X1991.02.18" "X1991.02.19" "X1991.02.20" "X1991.02.21" "X1991.02.22"
[421] "X1991.02.23" "X1991.02.24" "X1991.02.25" "X1991.02.26" "X1991.02.27" "X1991.02.28" "X1991.03.01"
[428] "X1991.03.02" "X1991.03.03" "X1991.03.04" "X1991.03.05" "X1991.03.06" "X1991.03.07" "X1991.03.08"
[435] "X1991.03.09" "X1991.03.10" "X1991.03.11" "X1991.03.12" "X1991.03.13" "X1991.03.14" "X1991.03.15"
[442] "X1991.03.16" "X1991.03.17" "X1991.03.18" "X1991.03.19" "X1991.03.20" "X1991.03.21" "X1991.03.22"
[449] "X1991.03.23" "X1991.03.24" "X1991.03.25" "X1991.03.26" "X1991.03.27" "X1991.03.28" "X1991.03.29"
[456] "X1991.03.30" "X1991.03.31" "X1991.04.01" "X1991.04.02" "X1991.04.03" "X1991.04.04" "X1991.04.05"
[463] "X1991.04.06" "X1991.04.07" "X1991.04.08" "X1991.04.09" "X1991.04.10" "X1991.04.11" "X1991.04.12"
[470] "X1991.04.13" "X1991.04.14" "X1991.04.15" "X1991.04.16" "X1991.04.17" "X1991.04.18" "X1991.04.19"
[477] "X1991.04.20" "X1991.04.21" "X1991.04.22" "X1991.04.23" "X1991.04.24" "X1991.04.25" "X1991.04.26"
[484] "X1991.04.27" "X1991.04.28" "X1991.04.29" "X1991.04.30" "X1991.05.01" "X1991.05.02" "X1991.05.03"
[491] "X1991.05.04" "X1991.05.05" "X1991.05.06" "X1991.05.07" "X1991.05.08" "X1991.05.09" "X1991.05.10"
[498] "X1991.05.11" "X1991.05.12" "X1991.05.13" "X1991.05.14" "X1991.05.15" "X1991.05.16" "X1991.05.17"
[505] "X1991.05.18" "X1991.05.19" "X1991.05.20" "X1991.05.21" "X1991.05.22" "X1991.05.23" "X1991.05.24"
[512] "X1991.05.25" "X1991.05.26" "X1991.05.27" "X1991.05.28" "X1991.05.29" "X1991.05.30" "X1991.05.31"
[519] "X1991.06.01" "X1991.06.02" "X1991.06.03" "X1991.06.04" "X1991.06.05" "X1991.06.06" "X1991.06.07"
[526] "X1991.06.08" "X1991.06.09" "X1991.06.10" "X1991.06.11" "X1991.06.12" "X1991.06.13" "X1991.06.14"
[533] "X1991.06.15" "X1991.06.16" "X1991.06.17" "X1991.06.18" "X1991.06.19" "X1991.06.20" "X1991.06.21"
[540] "X1991.06.22" "X1991.06.23" "X1991.06.24" "X1991.06.25" "X1991.06.26" "X1991.06.27" "X1991.06.28"
[547] "X1991.06.29" "X1991.06.30" "X1991.07.01" "X1991.07.02" "X1991.07.03" "X1991.07.04" "X1991.07.05"
[554] "X1991.07.06" "X1991.07.07" "X1991.07.08" "X1991.07.09" "X1991.07.10" "X1991.07.11" "X1991.07.12"
[561] "X1991.07.13" "X1991.07.14" "X1991.07.15" "X1991.07.16" "X1991.07.17" "X1991.07.18" "X1991.07.19"
[568] "X1991.07.20" "X1991.07.21" "X1991.07.22" "X1991.07.23" "X1991.07.24" "X1991.07.25" "X1991.07.26"
[575] "X1991.07.27" "X1991.07.28" "X1991.07.29" "X1991.07.30" "X1991.07.31" "X1991.08.01" "X1991.08.02"
[582] "X1991.08.03" "X1991.08.04" "X1991.08.05" "X1991.08.06" "X1991.08.07" "X1991.08.08" "X1991.08.09"
[589] "X1991.08.10" "X1991.08.11" "X1991.08.12" "X1991.08.13" "X1991.08.14" "X1991.08.15" "X1991.08.16"
[596] "X1991.08.17" "X1991.08.18" "X1991.08.19" "X1991.08.20" "X1991.08.21" "X1991.08.22" "X1991.08.23"
[603] "X1991.08.24" "X1991.08.25" "X1991.08.26" "X1991.08.27" "X1991.08.28" "X1991.08.29" "X1991.08.30"
[610] "X1991.08.31" "X1991.09.01" "X1991.09.02" "X1991.09.03" "X1991.09.04" "X1991.09.05" "X1991.09.06"
[617] "X1991.09.07" "X1991.09.08" "X1991.09.09" "X1991.09.10" "X1991.09.11" "X1991.09.12" "X1991.09.13"
[624] "X1991.09.14" "X1991.09.15" "X1991.09.16" "X1991.09.17" "X1991.09.18" "X1991.09.19" "X1991.09.20"
[631] "X1991.09.21" "X1991.09.22" "X1991.09.23" "X1991.09.24" "X1991.09.25" "X1991.09.26" "X1991.09.27"
[638] "X1991.09.28" "X1991.09.29" "X1991.09.30" "X1991.10.01" "X1991.10.02" "X1991.10.03" "X1991.10.04"
[645] "X1991.10.05" "X1991.10.06" "X1991.10.07" "X1991.10.08" "X1991.10.09" "X1991.10.10" "X1991.10.11"
[652] "X1991.10.12" "X1991.10.13" "X1991.10.14" "X1991.10.15" "X1991.10.16" "X1991.10.17" "X1991.10.18"
[659] "X1991.10.19" "X1991.10.20" "X1991.10.21" "X1991.10.22" "X1991.10.23" "X1991.10.24" "X1991.10.25"
[666] "X1991.10.26" "X1991.10.27" "X1991.10.28" "X1991.10.29" "X1991.10.30" "X1991.10.31" "X1991.11.01"
[673] "X1991.11.02" "X1991.11.03" "X1991.11.04" "X1991.11.05" "X1991.11.06" "X1991.11.07" "X1991.11.08"
[680] "X1991.11.09" "X1991.11.10" "X1991.11.11" "X1991.11.12" "X1991.11.13" "X1991.11.14" "X1991.11.15"
[687] "X1991.11.16" "X1991.11.17" "X1991.11.18" "X1991.11.19" "X1991.11.20" "X1991.11.21" "X1991.11.22"
[694] "X1991.11.23" "X1991.11.24" "X1991.11.25" "X1991.11.26" "X1991.11.27" "X1991.11.28" "X1991.11.29"
[701] "X1991.11.30" "X1991.12.01" "X1991.12.02" "X1991.12.03" "X1991.12.04" "X1991.12.05" "X1991.12.06"
[708] "X1991.12.07" "X1991.12.08" "X1991.12.09" "X1991.12.10" "X1991.12.11" "X1991.12.12" "X1991.12.13"
[715] "X1991.12.14" "X1991.12.15" "X1991.12.16" "X1991.12.17" "X1991.12.18" "X1991.12.19" "X1991.12.20"
[722] "X1991.12.21" "X1991.12.22" "X1991.12.23" "X1991.12.24" "X1991.12.25" "X1991.12.26" "X1991.12.27"
[729] "X1991.12.28" "X1991.12.29" "X1991.12.30" "X1991.12.31" "X1992.01.01" "X1992.01.02" "X1992.01.03"
[736] "X1992.01.04" "X1992.01.05" "X1992.01.06" "X1992.01.07" "X1992.01.08" "X1992.01.09" "X1992.01.10"
[743] "X1992.01.11" "X1992.01.12" "X1992.01.13" "X1992.01.14" "X1992.01.15" "X1992.01.16" "X1992.01.17"
[750] "X1992.01.18" "X1992.01.19" "X1992.01.20" "X1992.01.21" "X1992.01.22" "X1992.01.23" "X1992.01.24"
[757] "X1992.01.25" "X1992.01.26" "X1992.01.27" "X1992.01.28" "X1992.01.29" "X1992.01.30" "X1992.01.31"
[764] "X1992.02.01" "X1992.02.02" "X1992.02.03" "X1992.02.04" "X1992.02.05" "X1992.02.06" "X1992.02.07"
[771] "X1992.02.08" "X1992.02.09" "X1992.02.10" "X1992.02.11" "X1992.02.12" "X1992.02.13" "X1992.02.14"
[778] "X1992.02.15" "X1992.02.16" "X1992.02.17" "X1992.02.18" "X1992.02.19" "X1992.02.20" "X1992.02.21"
[785] "X1992.02.22" "X1992.02.23" "X1992.02.24" "X1992.02.25" "X1992.02.26" "X1992.02.27" "X1992.02.28"
[792] "X1992.02.29" "X1992.03.01" "X1992.03.02" "X1992.03.03" "X1992.03.04" "X1992.03.05" "X1992.03.06"
[799] "X1992.03.07" "X1992.03.08" "X1992.03.09" "X1992.03.10" "X1992.03.11" "X1992.03.12" "X1992.03.13"
[806] "X1992.03.14" "X1992.03.15" "X1992.03.16" "X1992.03.17" "X1992.03.18" "X1992.03.19" "X1992.03.20"
[813] "X1992.03.21" "X1992.03.22" "X1992.03.23" "X1992.03.24" "X1992.03.25" "X1992.03.26" "X1992.03.27"
[820] "X1992.03.28" "X1992.03.29" "X1992.03.30" "X1992.03.31" "X1992.04.01" "X1992.04.02" "X1992.04.03"
[827] "X1992.04.04" "X1992.04.05" "X1992.04.06" "X1992.04.07" "X1992.04.08" "X1992.04.09" "X1992.04.10"
[834] "X1992.04.11" "X1992.04.12" "X1992.04.13" "X1992.04.14" "X1992.04.15" "X1992.04.16" "X1992.04.17"
[841] "X1992.04.18" "X1992.04.19" "X1992.04.20" "X1992.04.21" "X1992.04.22" "X1992.04.23" "X1992.04.24"
[848] "X1992.04.25" "X1992.04.26" "X1992.04.27" "X1992.04.28" "X1992.04.29" "X1992.04.30" "X1992.05.01"
[855] "X1992.05.02" "X1992.05.03" "X1992.05.04" "X1992.05.05" "X1992.05.06" "X1992.05.07" "X1992.05.08"
[862] "X1992.05.09" "X1992.05.10" "X1992.05.11" "X1992.05.12" "X1992.05.13" "X1992.05.14" "X1992.05.15"
[869] "X1992.05.16" "X1992.05.17" "X1992.05.18" "X1992.05.19" "X1992.05.20" "X1992.05.21" "X1992.05.22"
[876] "X1992.05.23" "X1992.05.24" "X1992.05.25" "X1992.05.26" "X1992.05.27" "X1992.05.28" "X1992.05.29"
[883] "X1992.05.30" "X1992.05.31" "X1992.06.01" "X1992.06.02" "X1992.06.03" "X1992.06.04" "X1992.06.05"
[890] "X1992.06.06" "X1992.06.07" "X1992.06.08" "X1992.06.09" "X1992.06.10" "X1992.06.11" "X1992.06.12"
[897] "X1992.06.13" "X1992.06.14" "X1992.06.15" "X1992.06.16" "X1992.06.17" "X1992.06.18" "X1992.06.19"
[904] "X1992.06.20" "X1992.06.21" "X1992.06.22" "X1992.06.23" "X1992.06.24" "X1992.06.25" "X1992.06.26"
[911] "X1992.06.27" "X1992.06.28" "X1992.06.29" "X1992.06.30" "X1992.07.01" "X1992.07.02" "X1992.07.03"
[918] "X1992.07.04" "X1992.07.05" "X1992.07.06" "X1992.07.07" "X1992.07.08" "X1992.07.09" "X1992.07.10"
[925] "X1992.07.11" "X1992.07.12" "X1992.07.13" "X1992.07.14" "X1992.07.15" "X1992.07.16" "X1992.07.17"
[932] "X1992.07.18" "X1992.07.19" "X1992.07.20" "X1992.07.21" "X1992.07.22" "X1992.07.23" "X1992.07.24"
[939] "X1992.07.25" "X1992.07.26" "X1992.07.27" "X1992.07.28" "X1992.07.29" "X1992.07.30" "X1992.07.31"
[946] "X1992.08.01" "X1992.08.02" "X1992.08.03" "X1992.08.04" "X1992.08.05" "X1992.08.06" "X1992.08.07"
[953] "X1992.08.08" "X1992.08.09" "X1992.08.10" "X1992.08.11" "X1992.08.12" "X1992.08.13" "X1992.08.14"
[960] "X1992.08.15" "X1992.08.16" "X1992.08.17" "X1992.08.18" "X1992.08.19" "X1992.08.20" "X1992.08.21"
[967] "X1992.08.22" "X1992.08.23" "X1992.08.24" "X1992.08.25" "X1992.08.26" "X1992.08.27" "X1992.08.28"
[974] "X1992.08.29" "X1992.08.30" "X1992.08.31" "X1992.09.01" "X1992.09.02" "X1992.09.03" "X1992.09.04"
[981] "X1992.09.05" "X1992.09.06" "X1992.09.07" "X1992.09.08" "X1992.09.09" "X1992.09.10" "X1992.09.11"
[988] "X1992.09.12" "X1992.09.13" "X1992.09.14" "X1992.09.15" "X1992.09.16" "X1992.09.17" "X1992.09.18"
[995] "X1992.09.19" "X1992.09.20" "X1992.09.21" "X1992.09.22" "X1992.09.23" "X1992.09.24"
[ reached getOption("max.print") -- omitted 8498 entries ]
Wide to long data
library(lubridate)
#wide to long
SEUS_temp <- melt(allvalues_wide_latlon, id.vars = c("lat", "lon"))
colnames(SEUS_temp) <- c("lat", "lon", "date", "sst")
SEUS_temp[,date.edit := as_date(substr(date, 2, 11))][,year := year(date.edit)][,sst_avg := mean(sst, na.rm = T), by = year]
ggplot(data = SEUS_temp, aes(x=year, y = sst_avg)) +
geom_point() +
geom_smooth(method = lm) +
theme_classic()
Compare to SODA temp values
readRDS(here::here("Data","Spp_master","spp_master_ztemp.rds"))
readRDS(here::here("Data","Spp_master","spp_master_ztemp.rds"))